Possibility of "magic" trapping of three-level system for Rydberg blockade 
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The Rydberg blockade mechanism has shown noteworthy promise for scalable quantum compu- 
tation with neutral atoms. Both qubit states and gate-mediating Rydberg state belong to the same 
optically-trapped atom. The trapping fields, while being essential, induce detrimental decoherence. 
Here we theoretically demonstrate that this Stark-induced decoherence may be completely removed 
using powerful concepts of "magic" optical traps. We analyze "magic" trapping of a prototype 
three-level system: a Rydberg state along with two qubit states: hyperfine states attached to a 
J = 1/2 ground state. Our numerical results show that, while such a "magic" trap for alkali metals 
would require prohibitively large magnetic fields, the group IIIB metals such as Al are suitable 
candidates. 

PACS numbers: 32.10.Dk,32.80.Rm,37.10. Jk,31.15.A- 



Multiparticle quantum gates have been successfully 
implemented in neutral atoms using the Rydberg block- 
ade mechanism [ll-0| • In this method, qubits are encoded 
in hyperfine sublevels of the electronic ground state in 
each atom. Gate operations utilize interactions between 
pairs of neutral atoms with a separation on the order of 
10 iim. At this range, ground state atoms have negligi- 
ble interactions, while Rydberg states interact strongly 
with each other. Blockade refers to the inhibition of ex- 
citation of a target atom by the prior excitation of an- 
other (control) atom. This can be used to entangle the 
qubits and execute CNOT gates, which in principle forms 
a universal set of two-qubit gates for arbitrary quantum 
computation. 

Unfortunately, such a qubit implementation suffers 
from decoherence due to the trapping lasers. All three 
(two qubit and Rydberg) atomic levels are shifted via the 
Stark effect, leading to undesired differential phase accu- 
mulation. Even worse, the shift accumulation is uncon- 
trollable as it depends on the local laser intensity, which 
varies as the atom moves in the trap. 

This problem is similar to that encountered in optical 
lattice clock experiments. So-called "magic" traps were 
proposed as a powerful solution and have seen widespread 
use in clock experiments (see, e.g., reviews [^, Q and 
also 0, II)- Just as for neutral atom implementation 
of qubits, clocks depend on the transition frequency be- 
tween two particular atomic energy levels. With a par- 
ticular choice of trapping laser wavelength (and, in some 
cases, magnetic fields and polarization, see 0,11), both 
levels experience the same shift regardless of laser inten- 
sity, so the perturbation from the trap effectively van- 
ishes. Fig. [T] depicts the level structure. 

In the case of the alkalis, qubit transitions can be 
made Stark- and Zeeman-insensitive at a range of trap- 
ping wavelengths using a bias magnetic field tuned to a 
"magic" value, as described in The hyperfine struc- 
ture of Al is also favorable for "magic" trapping, dis- 
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FIG. 1: (Color online) Stark shifts affecting the Rydberg 
blockade three-level system. Qubit levels |0) and |1) are hy- 
perfine levels attached to the electronic ground state, while 
\R} is a Rydberg state. The unperturbed levels (solid) acquire 
a shift from the trapping laser, (a) For an arbitrary trapping 
laser, each state is shifted (dashed) by a different amount 
according to the local laser intensity, (b) With "magic" trap- 
ping conditions, all levels experience the same shift (dotted), 
and differential phase accumulation due to the trap vanishes. 



cussed in [9*1 and fllj, and the formalism from for 
alkalis is also applicable to Al. Thus, the problem of re- 
moving decoherence for the qubit levels has already been 
solved. 

However, excitation to a Rydberg level also causes 
decoherence of the qubit since in general the Rydberg 
trapping potential differs from the ground state poten- 
tial. This third level and additional source of deco- 
herence must be addressed. Previous proposals have 
considered "magic" trapping for a two-level system in- 
cluding a Rydberg state in alkalis [131 ^^.d in alkaline- 



earths [iSl, However, such treatments are not suf- 



ficient for the present problem, as the remaining Stark 
shifts of the qubit states would cause decoherence. We 
must therefore treat all three levels simultaneously. This 
problem is more akin to the situation in which con- 
sidered two-level systems in two species sharing the same 
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trap. 

Our goal in this paper is to find "magic" trapping con- 
ditions such that all three levels experience the same opti- 
cal trapping potential. We begin by reviewing the Stark 
effect theory underlying such "magic" traps. We then 
present results of numerical calculations for Rb, Cs, and 
Al. 

Non-magnetic states — We begin by summarizing the 
Stark effect formalism necessary for Rydberg states. This 
treatment is similar to that in 0, |l5j and references 
therein. The ac Stark shift of an atomic state |ni^, Mp) 



with total angular momentum F 
Mp can be written as 17 1 
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where £l and uj are the amplitude and frequency of the 
trapping laser field. The pre-factor a*°pM [oj) is the total 
polarizability of the state; note it depends on uj but not 
on £l. al^pMpi^) '^^Y be decomposed using irreducible 
tensor operators as 
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Here A is the degree of circular polarization (|^| < 1), 
while a^p, ctnF^ ^^'^ '^nF s-^^ the irreducible scalar, vec- 
tor, and tensor polarizabilities, respectively. The unit 
vectors are the laser wavevector (fc), laser polarization 
(e), andjDias magnetic field (B). Their relative geometry 
is as in [9|. The bias magnetic field is a static, externally 
applied field which defines the quantization axis. This 
"quantizing magnetic field" guarantees that Mp remains 
a "good" quantum number for the ac Stark effect pertur- 
bation formalism, from which Eq. ([2]) follows [3| • For lin- 
early polarized light, e ■ B — cos 6p, where 0p is the angle 
between the polarization and quantization unit vectors. 
If we consider circularly polarized light, defining e us- 
ing Jones calculus conventions, then |e • e^j^ = ^ sin^ 9k, 
where 9k is the angle between the wavevector and quan- 
tization unit vectors. 

Eqs. (HI and ^ determine our task: we must find uj 
and either 9k or 9p such that the Rydberg polarizability 
ctRydiuj) and the qubit state polarizabilities anFMpi^j) 
and OLnF'Mp, {uj) are all equal. Note that the polarizabil- 
ities in Eqs. ([T]) and ^ rnay be the conventional second- 
order quantities, as in [15j, or they may be replaced with 
third-order hyperfine mediated polarizabilities, denoted 
/3 and described in detail in The second-order def- 
inition considers only the interaction with the external 
electric field, neglecting the hyperfine interaction. This 
is sufficient for matching Rydberg and ground state po- 
larizabilities, but it cannot be used to find "magic" con- 
ditions for the qubit levels. As shown in Ref. [8|, if the 
hyperfine interaction is neglected, the qubit levels are 
degenerate and always experience the same shift, so any 
choice of lo and 9 is trivially "magic." Therefore, a third- 
order treatment including hyperfine and external electric 
field interactions is necessary for the qubit states. 

For Rydberg states, the tensor polarizability is highly 
suppressed at optical frequencies 12, [l^, so the total 



polarizability is dominated by the scalar part, which is 
essentially equal to that of a free electron, a^yci{oj) = 
— 1/cj^. Our ah initio numerical calculations confirmed 
this conclusion; see below for details. Since cx^y^ has 
negligible dependence on geometry (either 9k or 9p), the 
"magic" trap frequency will be determined by a^y^{Lij) — 
a*°,* (w), i.e., when the ground and Rydberg state po- 

1/2 

larizabilities are equal. 

Eqs. ([1]) and ([2]) also apply for qubit states in Al, with 
/3 in place of a as detailed above. Using Al in a linearly 
polarized trap, "magic" conditions for qubit transitions 
between nonmagnetic states will be set entirely by ge- 
ometry. For the alkalis, "magic" trapping is not possible 
using nonmagnetic states, so we consider magnetic states 
below. 

Magnetic states — Here we review the formalism to 
find "magic" conditions for qubit states with nonzero Mp 
projection. The methods developed above are not suffi- 
cient to find "magic" conditions for the qubit transition 
in alkalis. As was shown in [8], "magic" conditions do not 
exist for the nonmagnetic hyperfine states. This is due to 
the smallness of in comparison with a'^, along with 
the strict proportionality of a-nF'M , ^^'^ '^uFMf ■ How- 
ever, if we move to magnetic substates (with Mp ^ 0) 
we acquire unacceptable Zeeman sensitivity. 

A solution to this problem was developed in [lo| . 
Atoms with a J = 1/2 ground state are held in a cir- 
cularly polarized trap to take advantage of vector polar- 
izabilities. By utilizing multiphoton transitions between 
magnetic states with opposite projections (i.e., \nF' , Mp) 
and \nF, —Mp)) most of the first order Zeeman shift goes 
away as these states have opposite electronic g-factors 
(l9| . The remaining first-order shift is due only to the 
much smaller nuclear magnetic moment, which can be 
made to cancel the second-order shift with the applica- 
tion of a static magnetic field. The "magic" value of the 
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B-field is given by 



to second-order, it agrees well with the exact answer flo|. 



B„ 



2\{nF,MF,\fi-\nF',MF,)\- 



■ Wqubit , 



(3) 



where Wqubit is the energy splitting between the hyperfine 
levels. Although this expression for Bm is only accurate 
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Having dealt with Zeeman sensitivity, we turn to the 
Stark shift due to the trapping lasers. We need the differ- 
ential shift between the qubit states to vanish, i.e., both 
states must ex peri ence the same shift. This condition is 
satisfied when [10| 



WnF' -PnF) + + A COS Ok Mf 
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Here fS^ , /3" , and /3'^ are the third-order hyperfine me- 
diated polarizabilities referred to above, while ^ is 
the conventional second-order polarizability. Note the 
distinction between coupling schemes used in and 
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5/3'^ absorbs lengthy prefactors. The full form is 
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Bm is set by Eq. ([3|), while B is the actual applied field. 
Ideally B — Bm, but as suggested in this is not 
always possible. 

Clearly we need the third-order polarizabilities for the 
two states to cancel, but the appearance of the last term 
involving anSi/2 surprising. It arises due to inter- 

ference between the Zeeman shift and the vector part of 
the Stark shift, as these are both (axial) vector operators; 
see [l3| for a full discussion. Although /j.^ jib, a ^ 13, 
so this term is of a comparable order of magnitude and 
must be included in Eq. (|4]). 

Numerical evaluation — We used the same codes as 
in [§] for the ground states in all atoms considered. To 
summarize, we use the B-spline technique to generate a 
quasi-complete set of orbitals that are solutions to the 
Dirac-Hartree-Fock equations. To refine these solutions, 
we find the second-order self-energy operator to build the 
so-called Brueckner orbitals. Matrix elements are then 
calculated using the relativistic random-phase approxi- 
mation. 

For the Rydberg state calculations, we extended and 
modified our B-spline codes; the original codes are de- 
scribed in detail in [16[ . To generate a complete basis set 
including physically accurate Rydberg states, we dramat- 
ically increased the size of the cavity and the number of 
basis functions. As an illustration, a typical run to calcu- 
late low- lying states uses ~ 40 splines in a ~ SOas cavity. 
For calculations aimed at the 50s state, we obtained ac- 
curate results using ~ 200 splines in a ~ 8000a_B cavity. 
We also used a logarithmic rather than an exponential 
distribution of spline knots. This increased the accuracy 
of matrix elements by improving the representation of 



wavefunctions at large R near the cavity wall. 

Correlations were included by building the self-energy 
operator using a small basis set 40 splines in a 50aB 
cavity) and using this potential to build Brueckner or- 
bitals for a large set. This is justified since the self-energy 
operator diminishes rapidly outside the core, so highly 
excited states have a negligible contribution. Neglecting 
these states decreases calculation time dramatically. In- 
clusion of correlations introduced small but detectable 
corrections to Rydberg state energies and matrix ele- 
ments. Corrections were around the fourth significant 
figure for n = 50 states and diminished with increasing 
n. 

Results — We first present results for ®^Rb. This 
isotope has nuclear spin / — 3/2, and we are inter- 
ested in transitions between the |F' = 2,Mf' = 1) and 
\F = 1, Mf = —1) qubit states attached to the 5si/2 elec- 
tronic ground state. The hyperfine splitting is 6.83 GHz, 
and from Eq. ©, Bm ~ 3.25 G 0. In Fig. H we plot 
ctRyd, Q^fsjya' ^"^"^ ratio B/Bm- To achieve "magic" 
trapping using the scheme above, we require \B/Bm\ < 1 
in Eq. (g]). If \B/Bm\ is shghtly greater than 1, "nearly- 
magic" trapping is possible as considered in [10|. But in 
this case, \B/Bm\ diverges near Um, and the necessary 
B to achieve Stark-insensitive trapping is prohibitively 
large. The situation is qualitatively the same for ^■^'^Cs; 
our numerical calculations show a similar divergence near 
LOm- For the lighter alkalis, there are no doubly-magic or 
nearly-magic points for the qubit transition. 

Since it is not possible to build a "magic" three-level 
trap with the alkalis, we turn to ^^Al. As was shown 



ll| , "magic" trapping of its hyperfine states is aided 



by comparatively large vector and tensor polarizabili- 
ties. With nuclear spin / = 5/2, wc consider the \F' = 
3, Mfi = 1) and \F = 2, Mp — —1) qubit states attached 
to the 3pi/2 electronic ground state. The hyperfine split- 
ting is 1.506 GHz, and Eq. Q gives S,„ « 4.32 G. We 
may use either a linearly or circularly polarized trapping 
laser. In the case of a linearly polarized trap, the situa- 
tion is closely analogous to that in [9j . The "magic" wave- 



length is set where aRyd = ol 
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near w = 0.121 a.u. 
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FIG. 2: (Color online) Polarizabilities of the 5s state (dashed), 
the Rydberg state (solid), and the ratio B/Bm (lower frame) 
for the hyperfine transition \F = 2, Mf = 1) to \F — 1, Mf ~ 
— 1) in *^Rb. Since <C afs, the "magic" u simply occurs 
where a_R„d — af^ at approximately uj = 0.1062 a.u. (A — 
429 nm). B/Bm is obtained from Eq. Q. Near the circled 
"magic" Ld, B/Bm diverges, so magic trapping is impossible. 



(A — 377 nm). The enhanced tensor polarizability all 
but guarantees that the qubit transition can be made 
"magic." Numerical calculations show this does occur 
with a "magic" angle 9p ~ 65°. 

For a circularly polarized trap, no divergence of B/Bm 
occurs near the "magic" frequency LOm as for the alkalis. 
However, the "magic" condition is more complex, as the 
ground state second order a'^p^^^ is not negligible com- 
pared with afp^ ^. The total second-order polarizability 
of the ground state is given by Eq. ([2]) . This must 
equal the Rydberg polarizability, leading to 

aU^) + A cos0fe^a^(^) = a^M- (6) 

Our choice of uj and Ok must simultaneously satisfy 
Eqs. (HI and ([6]). We may solve for MpAcosOu appearing 
in Eqs. (j4| and ([6]) and equate the results, giving 
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FIG. 3: (Color online) Left-hand side (dashed) and right- 
hand side (solid) of Eq. (O. The range for "magic" trapping 
of the Rydberg transition is given by Eq. ([6]) and lies between 
the resonance at a; = 0.1155 a.u. (A = 394 nm) and the 
dotted line sX lo = 0.121 a.u. (A = 377 nm). The "magic" 
frequency and "magic" angle for the qubit transition are set 
by the curves intersection, just under a; = 0.121 a.u. This 
combination of iOmagic and Omagic will allow Stark and Zeeman 
insensitive trapping for the three-level system in Al. 



We plot the left and right hand sides of Eq. ([7]) in 
Fig. [3l Since they intersect in the range allowed by 
Eq. © , "magic" trapping for this three-level system in Al 
is possible. While circularly polarized trapping would be 
more complex than the linearly polarized trap presented 
above, it has the additional advantage of less Zeeman 
sensitivity. This is because the "magic" i3-field removes 
Zeeman effects to second-order, while the linearly polar- 
ized trap only removes Zeeman decoherence to first order 
in the B-field. 

We have presented a method to remove differential 
Stark shifts for a three-level atomic system consisting of 
a Rydberg state and two hyperfine states attached to 
the ground electronic state. Such system is an essential 
element of the CNOT gates utilizing Rydberg blockade 
mechanism. Although such a trap is not possible for 
the alkalis, our numerical calculations show that Al may 
be trapped using this method. Such a trap could prove 
useful for removing decoherence from trapping lasers in 
implementing the Rydberg blockade mechanism. 
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